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I. INTRODUCTION 

Since the recent discovery of the strongly correlated 
QCD fluid, the phenomenology of quark-gluon plasma 
(QGP) has attracted a lot of attention [1]. Many pre- 
dictions have been made for the several possible phases 
of the QGP in which quarks and gluons arc no longer 
confined to nucleons and pions [2-5]. At extremely high 
temperature and low density, thermal excitations of the 
gluon field are expected to screen interactions between 
color charges and ultimately, due to asymptotic freedom 
result in a weakly interacting quark gas [6-8]. Simi- 
larly, at low temperature but high density Debye screen- 
ing is expected to reduce the range of strong interac- 
tions. Novel phases that correlate quark color and fla- 
vor are predicted to occur at asymptotic densities due to 
the attractive nature of quark-quark interactions in cer- 
tain color-flavor locked combinations [9, 10]. Thus finite 
temperature and/or density are expected to reflect on 
various aspects of confinement. Furthermore any mod- 
ification of the ground state influences symmetry prop- 
erties and in particular restoration of chiral symmetry 
is expected. At finite density and low temperature the 
precise relation between chiral symmetry restoration and 
dcconfincmcnt is not yet known. A common wisdom 
is that in this regime there is a phase transition from 
hadronic (confined) matter at low density to the uncon- 
fined (possibly superconducting) phase at higher densi- 
ties [11-14]. Recently, however it has been observed that 
the confined and deconfmed phases may be separated 
by a phase where quarks are confined but chiral sym- 
metry is restored. In the limit of an infinite number of 
colors Nc — ► oo this so called quarkyonic phase would 
in fact extend to infinite density, since in this limit De- 
bye screening due to quark loop vanishes [15]. In this 
paper we examine the possible emergence of this new 
phase using a canonical formulation of the QCD many 
body problem in the Coulomb gauge. The Coulomb 
gauge canonical formulation can describe both finite tem- 
perature and density. In this formulation manifestation 
of deconfmement can be inferred from the temperature 
and/or density dependence of the color Coulomb inter- 
action. In contrast lattice simulations at finite density 
arc still at their infancy [16, 17]. One of the smoking 
gun signals of dcconfincmcnt is the possible dissociation 



of heavy quarkonia [18-25]. By studying density depen- 
dence of the Coulomb gauge heavy quark potential we 
will be able to explore charmonium properties at finite 
density. 

The paper is organized at follows. In the following 
section we discuss the Coulomb gauge QCD and approx- 
imations relevant to the problem in hand. Zero and fi- 
nite density properties are discussed in Sections III and 
IV, respectively Charmonium dissociation is studied in 
Sec. V and followed by conclusions and outlook. 



II. COULOMB GAUGE QCD 

In this section wc briefly discuss QCD in the Coulomb 
gauge and the approximations appropriate for the high 
density and/or temperature systems [26, 27]. In the 
Coulomb gauge gluons are described by the transverse 
potentials, A a (x), a = 1 • • • N% - 1, V • A a (x) =0 and 
the conjugated, transverse momenta, II a (x), 

[A a (x),Il b (y)]=zS T (x-y)S ab , 



(1) 

where <5 T (x - y) = [I - W/V 2 ](5 3 (x - y). The canon- 
ical momentum II a (x) is the negative of the transverse 
component of the chromo-electric field. The quark and 
antiquark degrees of freedom will be defined below in 
terms of the canonical set of Dirac fields, ^j(x), ^j(x), 
i = 1---Nc, one for each flavor. The Hamiltonian is 
given by 



H = Hd + Hym + He 



(2) 



where Hp contains the Dirac kinetic energy and quark- 
transverse gluon interaction, Hym is the Yang-Mills 
term, which contains the gluon kinetic energy and the 
three- and four-gluon interactions and finally He is the 
Coulomb potential given by, 

He = \J d^dyJ- 1 p a ( X )K(^a;y,b)[A}Jp b (y). (3) 

It represents the non-abelian Coulomb gauge interaction 
between color charge densities, /O a (x) = ^ {x)T a ij){x) + 
/ a t c A b (x)n c (x), mediated by the Coulomb kernel K[A] 
given by 



K(x,a;y,b)[A] = 



V -V 



(-v 2 



(x,a;y,b) 



(4) 
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Finally J = Det [—V • 27] is the determinant of the 
Faddeev-Popov operator; V = V ab = <5 a &V + <?/ ac bA c is 
the covariant derivative in the adjoint representation, and 
B is the chromo- magnetic field, B a (x) = V x A a (x) + 
(fl/2)/nkA i (x)xA c (x). At every space point-x and color 
component a the Coulomb gauge potentials, A a (x) are 
an analog of a curvilinear coordinate. This is because 
their values are restricted to reside within the bound- 
ary of the Gribov region which has a nontrivial metric 
determined by J [28]. A confinement scenario in the 
Coulomb gauge states that it is the field configurations 
near the boundary of the Gribov region <9f2 that domi- 
nate the QCD vacuum [29, 30] and it follows that fluc- 
tuations near the boundary lead to massive quasiparticle 
excitations. Since the exact parametrization of the Gri- 
bov horizon is not known the quantitative description of 
this confinement scenario varies depending on how re- 
striction to the Gribov horizon is implemented, but the 
general features seem to be robust [31-34]. With this pic- 
ture in mind we approximate the ground state of a finite 
density quark plasma by the state with no quasiparticle 
gluon excitations and the large background fields con- 
centrated near dfl lead to enhancement in the long range 
behavior of the Coulomb kernel [26, 34]. That is we make 
the replacement 

iT(x,a;y,&)[A] -► (K (x, o;y, 6) [A]} =K(x-y)5 ab (5) 

with the potential K(r) modified from its free (A = 0) 
form K{r) = a/r, for large r due to the large fields A £ 
dfl. In particular we approximate the kernel by the from 



K{r) = K c (r) + K L {r) 



(6) 



with Kq and Kl being the short-range Coulomb and 
long-range linear potentials, respectively which will be 
discussed in detail in the following section. The final 
Hamiltonian, describing massless quarks with energies 
below gluon quasiparticle excitations is given by 



ff = /^t (x)Ha . VM ,) 



+ - J dxdy^(x)T a ^(x)K(\x-y\)^(y)T a ^(y). 

(7) 

Note, that from Eq. (7) it follows that it is —CfK(t) 
which is the instantaneous interaction in the color-singlet 
qq channel. 



III. QUARKS AT ZERO DENSITY 

In presence of the effective density-density interaction 
mediated by the kernel K, in Eq. (7) quarks and an- 
tiquarks acquire effective mass, which in the mean filed 
approximation can be described within the Hartree-Fock- 
Bogolubov framework [11-13, 35]. The single quark 



quasiparticle operators are defined by a canonical trans- 
formation to a plane wave representation of the Dirac 
fields. 



E 

A 



J dke i;k x [ U (k, A)6 k . A + w(-k, A)4 >A ] 



(8) 



and similarly for Here dk = d 3 k/(2n) 3 : A = ±1/2 
is the quark (antiquark) spin projection and b(d),w(d*) 
are the quark (antiquark) annihilation and creation oper- 
ators, respectively. These quasiparticle operators satisfy 
the standard fermion anti-commutation relations and de- 
fine the vacuum state, by &A,k|vac) = d\^\va,c) = 0. The 
single particle wave functions are given by 



u T (A,k) 



v 1 (A,k) = -= 




o- • kx/3 



' 1 - — f ' kX/9, \/l + — Xf3 
m k V rn k 



(9) 



where x — *°2X) m k/E k = sin^ and is the BCS 
angle, which determines the number density of quark- 
antiquark pairs in the BCS vacuum, i.e. the quark con- 
densate and the extent of chiral symmetry breaking [36- 
38]. At zero temperature and density the BCS angle is 
determined by minimizing the vacuum energy density, 
5(vac|JT | vac) = 0, which leads to the gap equation 



sin < 



COS( 







P*p = % / dkK(\k 



p\)[skC p - c k s p k ■ p] 



(10) 



Here K(p) is the Fourier transform of the effective po- 
tential from Eq. (4). The Fourier transform of the linear 
potential, Kl has to be taken with care, since naively, 
/ dx|a;| exp(— k ■ x) = oo. We introduce an infrared reg- 
ulator, e and define [39] 



* lW ^, (r ,=!(i-^ 



e 



so that lim c _»o K^ ^r) = — br. The difference between 
the liner potential and the IR finite approximation is 
shown in Fig. 1 . In momentum space the IR finite kernel 
becomes 



KlAp) 



8-Kb 



p 2 (p 2 



-^(2.) 3 <5 3 (p) (12) 



It is clear that the S term does not contribute the gap 
equation (10) and the gap equation is well defined in the 
limit e — > 0. Equivalcntly, the gap equation is invariant 
under a constant shift in the potential, 



K{r) -> K(r) + C. 



(13) 
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Since such a shift induces a contribution to the Hamilto- 
nian proportional to the square of the total charge oper- 
ator J2 a Q a Q a , 



Q a = J dxp a (x) 



(14) 



invariance under (13) is an exclusive property of color sin- 
glet states and it is only matrix element invariant under 
global color rotations that are are physical. In contrast, 
under this shift, the energy f2 p of a single quark state ( 
&t(A, p)|vac)) that is given by, 



n. 



pc p 



2 



dkfif(|k- p\)[s k s p + c k Cpk ■ p] (15) 



transforms to 



C F 

flp > flp -\ — C. 



(16) 



So single quark states are clearly unphysical. For e — > 
the 5-term in K dominates the integrand in Eq. (15) and 
the quark self-energy becomes negative and tends to — oo 
in the e — > limit. As pointed out in [39] this is necessary 
in order for color-singlet qq excitations to have finite, non- 
negative energies. This is because the potential energy in 
the qq bound state, given by ~CfK(p), where p is the rel- 
ative momentum between the quark and the antiquark is 
large and positive for small e (and approached +00 in the 
limit e — > 0). Thus the infinities in the e — > limit can- 
cel between the self-energies and the residual interaction 
between the quark and the antiquark. Since (f> 7^ is a 
lower energy state compared to <f> = 0, after cancelation 
of the IR divergencies, the finite energy of color single 
excitations is C-indcpendent and non-negative. This is, 
however not the case for color-nonsinglet states. For ex- 
ample a single quark state with energy given by Eq. (15) 
has negative energy which becomes —00 in the limit when 
the interaction is confining (i.e. e — > 0). This is clearly 
unphysical as one would expect colored states to have 
positive, and IR diverging energies, Q — > +00 in the con- 
fining limit. Since the shift in Eq. (13) is a symmetry 
of the physical sector we can choose C to cancel the IR 
divergence in the quark self energy, and redefine Kl ac- 
cordingly Thus instead of Eq. (12) we should use 



Kl.M 



8tt& 



p 2 (p 2 



(17) 



which is positive, gives f2 — > +00 as e — > and does not 
affect properties of color singlet states. In fact and inter- 
action without the IR divergent constant term constant 
term is obtained when computing the expectation value 
of K[A] (cf. Eq. (4)) in a mean filed ansatz for the gluon 
vacuum distribution [26, 27]. 

A similar argument for regulating the p = singular- 
ity of Kl was proposed in [40]. There by explicitly re- 
stricting the spectrum of the Hamiltonian to include only 
the color-singlet subspace an even stronger constraint on 



the IR momentum dependence of the kernel was derived, 
namely Kp(0) = 0. With an IR regulated (finite-e) kernel 
the minimal physical requirement, however is that col- 
ored excitations have positive energies, (becoming +00 
as e — > 0). Thus we do consider kernels K(p) which have 
intcgrablc, but in principle finite zero-modes. 



V(r) = b r 




FIG. 1: IR finite approximation to the linear potential of 
Eq. 11. The proper potential for describing color-nonsinglet 
state, (cf. Eq. (11)) corresponds to a downward-shift by two 
units. 



IV. QUARKS AT FINITE DENSITY 

With the interactions in the Hamiltonian now well de- 
fined in the IR, the system at finite quark density can be 
described using standard many body techniques [11-13]. 
In particular the finite density gap equation becomes, 



psp 



2 



d\iF k K eff (\\t- p|)[s fe c p - c fc Spk ■ p] (18) 



and the single quark energy is given by 
C F 



~PC P + — I dkF k K eff (\k - p\)[s k s p + c fc c p k • p]. 



(19) 

Here F k = l—n k ~n k and n k (n k ) is the quark (antiquark) 
occupation number at zero-temperature, 



n k = 0(0* -k F ),n k = 6(Q k + k F ) 



(20) 



with kp denoting the the Fermi momentum. Because of 
the Pauli blocking factor F k interactions between quarks 
have small effect on quark levels inside the Fermi sphere 
where quarks arc effectively free f2 p ~ pc p . The self- 
energy contributes mainly for states above the Fermi sur- 
face and results in states which are similar to the confined 
quark states at zero-density. Since any modification of 
the antiquark distribution from its vacuum value is sup- 
pressed at finite quark density, it will be ignored in the 
following and we will only consider quark-hole excitations 
near the quark Fermi surface. 
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A. Screening effect at finite density 

At finite density the quark-quark potential CfK(p) 
is screened by particle-hole excitations near the Fermi 
surface, resulting in the effective interaction CpK e f f(p) 
which enters Eqs. (18), (19), 



Keffip) 



K(q) 
K{q) 



Keff{q) = 



K(q)IL(q)K(q) + ■ ■ 
K(q)H(q)K eff (q), 

K{q) 



1 + K(q)Tl(q) 



(21) 



(22) 



The vacuum polarization, 11(g) is shown in Fig. 2 and it 
describes the probability for creating a particle-hole ex- 
citation. Since the pair is excited in the colored state in 
the confining limit (e — > 0) it is expected that H(q) — > 0, 
since the amplitude is proportional the the inverse of the 
particle-hole excitation energy, which becomes infinite in 
the confining limit. Assuming phase transition does take 
place, however we can examine the effect of Dcbye screen- 
ing in the deconfined phase near or above the phase tran- 
sition density. In this case vacuum polarization is given 
by 



n(«) 



2 



n k - n| k+q | 



r> fe -n 



|k+q| 



(23) 



1 + SfcSjk+ql + CfcC| k+q 



k • (k + q) 



fc(|k + q|) 



where nf is the number of light flavors. The set of equa- 
tions, (18), (19), (22) and (23) with K(p) given by 



k( P ) 



4ira(p) 



8nb 



p2 {pi _|_ £ 2 ^ 



(24) 



forms a set of coupled equations which we solve numeri- 
cally and discuss below. 



quarks flavors using 
a(p) = 



AirZ 



QCD 



(25) 



with Z = 5.94, c = 40.68 and A QC d = 250 MeV deter- 
mined from fitting the zero-density qq potential [41]. At 
zero density what makes quark energies finite is the IR 
regulator e in Eq. (24), at finite density, above the decon- 
finement phase transition we expect, however, that the 
self-consistent set of equations will admit nontrivial so- 
lutions in the limit e — > 0, with the Fermi momentum 
taking over the role of the IR regulator instead. We have 
verified numerically that this indeed is the case. The 
equations are solved by iterations. For given density, k F 
we start with a small, but finite e, e.g. e = 0.008 GeV, so 
the self-energy at the first iteration has a sharp jump at 
the Fermi surface (see left panel in Fig. 5). At this ini- 
tial state the vacuum polarization is highly suppressed 
(see left panel in Fig. 6), particularly for small momenta. 
After a few iterations, however, the self-energy becomes 
regulated by the vacuum polarization itself and simulta- 
neously the vacuum polarization increases at small mo- 
menta. Finally, we reduce the initial value of the IR 
regulator, e and we verify that after several iterations 
solutions converge to the same value, regardless of the 
starting value of the regulator. We repeat the calcula- 
tions for several values of the Fermi energy. We show 
these finial results for the effective potential, BCS gap 
angle, quark single energy and vacuum polarization in 
Figs. 3, 4, 5, 6, respectively. 

It is often stated that chiral symmetry restoration and 
deconfinement occur simultaneously [11, 12]. Our cal- 
culation illustrates that this need not be the case. At 
finite density, effective potential is already deconfined, 
but the gap equation admits nontrivial solutions. Only 
when the Fermi momentum increases above, approxi- 
mately kp ~ 0.05 to 0.06 GeV, the effective potential 
at large distance is not strong enough to sustain sponta- 
neous chiral symmetry breaking. 



3. 



anti-Q 

FIG. 2: Vacuum polarization fl 9 . 



B. Numerical result 

We solved coupled equations (18), (19), (22) and (23) 
at finite density and zero temperature with two light 



C. Quarkyonic matter 

Another manifestation of the independence of chiral 
symmetry restoration and deconfinement is the appear- 
ance of confined, chiraly symmetric, quakyonic, mat- 
ter [15]. Since K oc g 2 oc Nc we can extract explicit 
Nc dependence of the effective interaction, 



_3_f>N c =3 
j^-Nc _ N C 



(26) 



The qq potential being proportional to CpK^y at large 
Nc becomes 



C F K eff _ \k»°=* 



(27) 
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FIG. 3: Effective potential at finite density (ro 
l/OASGeV- 1 ). 



FIG. 5: Left panel : Variation of single quark energy Q p 
through iterations, for k,F = 0.3 GeV. f2 p is the largest at 
first iteration. Right panel: Q p at the end of the iterations. 





FIG. 4: Solution of gap equation <f>k at finite density. For 
k,F > 0.05 GeV the solution of the gap equation is 4>{p) = 0. 



FIG. 6: Same as in Fig. 5 for vacuum polarization II 9 . In the 
left panel, Tl(q) is the smallest at first iteration. 



since the vacuum polarization contribution is suppressed 
in the Nc — > oo limit. In this limit we thus find that the 
Debye screening disappears and confinement is restored 
at any density. The gap equation (18), however is not 
affected by the large Nc limit, and the Pauli blocking 
remains in effect. Thus as density increases it will even- 
tually prevent gap equation from developing nontrivial, 
charily broken solution [14]. 



CHARMONIUM BINDING AT FINITE 
DENSITY 



Modifications of charmonium properties, in particular 
binding energy and size may be strongly affected when 
the bound state propagates through the plasma. At fi- 
nite temperature such modifications have been studied by 
computing charmonium spectrum using a temperature- 
dependent static potential that is extracted from the lat- 
tice calculation [18-20]. With such temperature depen- 
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dent static potential, melting of charmonium can be in- 
ferred from the temperature dependence of the spectral 
function [24, 25] or by directly solving the bound state 
Shrodinger equation [22, 42]. In the former approach, 
melting of charmonium can be seen when the bound state 
peak of spectral function collapses and as temperature 
increases it becomes buried under the continuous back- 
ground. From Shrodinger equation calculation melting 
is inferred at a temperature when the bound state solu- 
tion disappears. Melting temperature of the J/vP, from 
these two different approaches appears to be in a very 
good agreement at approximately 1.6T C . Here we use the 
Shrodinger equation approach to calculate dissociation of 
the J /if? at finite density with two flavors of light quarks. 
The mass Mjpc of the J /if? is given by 

Mjpc = ejpc + 2m + C F [K eff {r ) - K CL (r )], (28) 

and 

(ejro-£)*Z{q) (29) 
m 

= -^E / rfq'^(q-q')^ e //(|q-q'l)^(9')- 

a' 

with N being the radial quantum number, a = (S q ,L q ) 
stands for total quark spin S q and relative orbital angular 
momentum L q . The heavy (charm) quark mass, m ab- 
sorbs any finite shift that has been removed from K e f f as 
discussed in Sec. Ill and is fixed, by fitting spin-averaged 
L q = charmonium masses M^ s = i[M -+ + 3Mi— ] = 
3.068 GcV. The energy of J /if? and its first radial exci- 
tation as a function of the light quark density are given 
in Table I and the binding energy is plotted in Fig. 7. 
From our numerical evaluation is follows that melting 
occurs at light quark density kp ~ 0.9 GcV for J /if? and 
kp ~ 0.3 GeV for Finally in Fig. 8, we show how the 
wave function of J /if? collapses with increasing density. 



TABLE I: J/if! and energy spectrum at finite density 



fc F (GeV) 


J/*(GeV) 


*'(GeV) 


< r > Jf y (GeV- 1 ) 


0.0 


3.065 


3.823 


1/0.476 


0.2 


2.801 


3.288 


1/0.417 


0.4 


2.771 


2.946 


1/0.351 


0.6 


2.733 


2.771 


1/0.250 


0.8 


2.682 


2.683 


1/0.050 



VI. CONCLUSION AND OUTLOOK 

In this work we studied the screened quark-quark ef- 
fective potential at finite temperature and presented nu- 
merical result for the effective potential. We investigated 
the restoration of chiral symmetry using the many body 
framework of the Coulomb gauge QCD and find that 




,i i i i i i i i i i i 

0.2 0.4 0.6 0.8 1 

k F (CeV) 



FIG. 7: ej/q, and e*' at finite density, where Ejpc is defined 
in Eq.(29). 




r(GeV ') 



FIG. 8: J/*'s Wave function at finite density. 



the transition densities of chiral symmetry restoration 
and de-confinement need not be related. In the decon- 
fined phase, chiral symmetry can be broken at density 
kp ~ 0.05 to 0.06 GcV. We also computed the melt- 
ing density for J /if? and if?', which we find to be at 
kp ~ 0.9 GcV and kp ~ 0.3 GeV respectively. Simi- 
lar to the situation at high temperature and low density 
[21, 22, 43, 44], collision of gluons and quarks with heavy 
quarkonium will reduce the dissociation temperature of 
quarkonium. We are expecting the same situation oc- 
curs at high density and low temperature. As pointed 
out previously [21, 22, 43, 44], production of quarkonium 
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in heavy-ion collision may provide clear signal for quark- this case, the detail discussion of collision dissociation of 
gluon plasma. Clearly, the collision effect of gluons and heavy quarkonium at the phase of high density and low 
quarks with heavy quarkonium should be incorporated in temperature is under way. 
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